Phase slip phenomena in superconductors: from ordered to chaotic dynamics 
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We consider flux penetration to a 2D superconducting cylinder. We show that in the low field limit 
the kinetics is deterministic. In the strong field limit the dynamics becomes stochastic. Surprisingly 
the inhomogeneity in the cylinder reduces the level of stochasticity because of the predominance of 
Kelvin-Helmholtz vortices. 
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The kinetics of vortex production in superconductors 
and superfluids is one of the intriguing problems of con- 
densed matter physics. It is interesting not only in the 
field of solid state physics but represents as well a very 
good model to study topological phase transitions in cos- 
mology and other branches of physics [l[ . In the last few 
decades different scenarios for vortex production in su- 
perconductors and superfluids were proposed. The most 
common way to produce vortices is to increase the su- 
pcrfluid velocity in order to reduce the energy barrier 
between homogeneous flow and flow with vortices. This 
mechanism is observed in rotating 3 He where vortex nu- 
cleation and critical velocities are measured [3j. In 2D 
homogeneous superconducting films, increasing the cur- 
rent leads to dynamics which are similar to the phase 
slip (PS) transition in ID [4]. The order parameter (OP) 
reaches zero along a straight line across the film and the 
phase displays a 2tt jump along this line. This PS line 
solution |5[ corresponds to the deterministic and most 
ordered PS kinetics in 2D. The inhomogeneity caused by 
current contacts leads to a qualitatively similar picture. 
The OP is strongly suppressed along a straight line across 
the film but it reaches zero only at two points on this line. 
This pair of vortices is called kinematic vortex-antivortex 
(VaV) pair [f|. It spreads quickly in opposite directions 
along this line propagating the 27r jump of the phase. 
Therefore PS occurs without formation of well defined 
VaV pairs Q. 

A different scenario of vortex production was proposed 
by Kibble[|] and Zurek@ (KZ). When the sample is 
quickly quenched through the critical temperature T c , 
the nucleation of the low temperature phase starts in 
different places with uncorrelated phases of OP. Then, 
domains grow and start to overlap leading to the forma- 
tion of vortices. This mechanism is a promising way to 
test cosmological theories in condensed matter physics 
[ToL 11 1 . This dynamics is stochastic and sensitive to 
small variations of initial conditions. On the other hand 
the dependence of the vortex density on the quench time 
and their spatial correlation are universal. Later, in Refs. 
12l 13 1, it was proposed that the quench occurs not only 
due to fast temperature change but also due to the tem- 
perature front propagation. Aranson et al. considered 
the case of a temperature quench in the presence of ex- 




FIG. 1: Geometry of the system : a 2D cylinder with an 
applied magnetic field H. 



tcrnal current 14J. The new phase with zero current 
grows after the quench. Therefore on the border of the 
quenched region, the superfluid velocity has tangential 
discontinuity, leading to vortex formation, similarly to 
the classical hydrodynamic Kelvin-Helmholtz (KH) in- 
stability [1_5| which is also known in superfluids 
Moreover, the KH instability suppresses the development 
of KZ vortices [14j . In this paper, kinematic VaV, KZ 
and KH vortices are distinguished by their production 
mechanism although they are topologically equivalent. 

To demonstrate how deterministic type of dynamics 
becomes stochastic, we model a superconducting film 
rolled on a cylinder in an external time dependent mag- 
netic field parallel to the cylinder axis (Fig.l). Depending 
on the applied magnetic field and the dimensions of the 
ring, we follow the evolution from the deterministic PS 
line dynamics to the stochastic behavior described by the 
KZ mechanism. In the proposed model, topological de- 
fects arc generated by the intrinsic quench induced by 
the external field. The evolution towards stochastic be- 
havior is strongly influenced by the KH instability which 
develops in the presence of inhomogeneities. To model 
the inhomogeneity of the film we assume that there is a 
thin stripe of superconductor along the film with a differ- 
ent coherence length. The thickness of the film d is small 
d <C £ -C A c ff. Here £ is the coherence length and A e ff 
is the Pearl penetration depth. Therefore we can neglect 
all corrections to the external magnetic field H caused by 
the current in the film. The radius of the film is R > £. 
The time dependent Ginzburg-Landau (TDGL) equation 
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in dimensionless units has the form: 

+ = K*)^ " V#| 2 ) ~ (*V + a) 2 V + 77. (1) 

Here ip is the dimensionless complex OP, the spacial co- 
ordinate r is measured in units of £ and time is measured 
in units of phase relaxation time Tg = 47rA e f f(T " , X e s = 4r , 
A is the bulk penetration depth, a n is the normal state 
conductivity, and c is the speed of light. The parameter 
u = ^ is a material dependent parameter, where Ty, is 
the relaxation time of the amplitude of the OP. According 
to the microscopic theory, u is ranging from 5 to 12 but 
we assume < u < 00. The vector potential a is mea- 
sured in units of ^7 where 4>o is the flux quantum. The 
function b(z) models the z dependence of the coherence 
length £(z). As shown in Fig.l, we chose £(z) = £/iJb(z) 
and b(z) = 1 — b 2r d(z - iw/4)i9(3it;/4 - z). Here w is the 
width of the film in units of f , b parameterizes the level 
of inhomogeneity of the film and -d(x) is the Heaviside 
step function. Here we use periodic boundary conditions 
and the boundary condition with vacuum [18j at z = 
and z = w. The equation for the electrostatic potential 



<I> measured in units of 



2TTCT0 



, where e is the electronic 



charge and H is the Planck constant, reads: 
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To model the process of vortex formation we assume that 
at time t < external magnetic field is absent. At t = 
the field suddenly appears and stays constant for t > 
i.e. tangential component of the vector potential is a/S{t). 
We thus study the kinetics of the vortex generation as a 
function of a with different values of u. 

Let us first consider the stability of the solution in the 
uniform case. We linearize the TDGL Eqs. (I1I2[ ) in small 
fluctuations of OP /(r, t) = ip(r,t) — ipo and search for 
a solution in the form f(r,t) — Ck exp (ikr + X^t). 
It is clear that the transverse k z component always con- 
tributes to the stability of the initial state. Therefore, the 
condition Ak > is the same as in ID » 
where (f> is the magnetic flux through the ring at t > 0. 
This condition provides a rough estimate for the number 
of the expected PS events N ~ 4^. It defines the first 

critical value of the external field a c i = 1/V3- There- 
fore, in the low field limit a c \ < a, < 1, the dynamics will 
be similar to the ID case with very weak z-dependence. 
Well defined vortices may appear in this region of the 
field if the film is inhomogeneous as shown in Fig.l. The 
situation is different when the field a increases further. 
Dropping k z — 0, the eigenvalues are: 
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describes the decay rate of the uniform solution. On 



the other hand for finite k, Ak is positive and char- 
acterizes the growth of the corresponding Fourier com- 
ponents Ck- The fastest growth is found for k = 
—16utp 2 a 2 — tPq(u — 2) 2 + 16a 4 and the rate is de- 
termined by A max = Ygi^r (8(u - 4)^oa 2 + 16a 2 + ^(u- 
2) 2 ). The qualitative difference in kinetics takes place 
when the decay rate of the uniform solution becomes 
faster then the growth of the new phase. This effect is 
similar to the quench through T c in the KZ mechanism 
[J 13]. We find that at a > a c2 = V2, the OP is sup- 
pressed to zero and the growth of the phase with finite 
k is accompanied by the rapid development of vortices. 
The density of vortices may be estimated using Zurck 
arguments where the quench time should be replaced by 
tq = (a 2 — 1) _1 leading to n cx Tq 1 ^ 2 [li|. 
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We simulate Eqs. (|ll2[ ) using the fourth order Runge- 
Kutta method. The spatial derivatives are evaluated us- 
ing a finite difference scheme of second order or using a 
fast fourier transform algorithm depending on the bound- 
ary conditions. The choice of the algorithm is made 
to optimize the convergence and the calculation times. 
The calculations are performed for the vector potential 
< a < 5 and for the total flux tfi through the ring 
ranging from to 50<fo- 

We investigate the flux penetration into the homoge- 
neous ring for two different boundary conditions. In the 
case of periodic boundary conditions we identify different 
regimes in accordance with Eq. (|3]) . In the small field limit 
a < a c i the ring is in a stable state and the penetration 
of the magnetic flux into the ring can only be induced 
by a very strong noise rj in Eq.([T]). When a c \ < a < a C 2, 
in agreement with the stability analysis, the PS kinetics 
depends on the external magnetic field. When ^ < 10 
the kinetics is similar to the ID case. The transition 
is characterized by one or more lines in the z-direction 
where the OP decreases to zero (the PS line case |6J]). 
These lines may appear simultaneously or consecutively 
in time, depending on u 4]. As expected, the number 
of PS events is determined by the ratio These PS 
lines represent the limiting case of kinematic VaV pairs 
travelling with infinite velocity. 

When the flux is increased (4- > 10), the kinematic 
vortices become clearly distinguishable. In Fig. 2(a), we 
present the time evolution of the average value of the 
OP together with the time dependence of the number 
of vortices in the sample. The kinetics is characterized 
by series of consecutive PS events well separated in time 
(Fig. 2(a)). As it was noticed [23], few VaV pairs may 
propagate along the same line at the same time. PS 
events are produced by kinematic VaV pairs propagat- 
ing along the same line where the amplitude of OP is 
reduced. Kinematic vortices can propagate in the same 
direction, one after another or in opposite direction lead- 
ing to annihilation of VaV pairs and accelerating the dy- 
namics. Contrary to 0, kinematic VaV pairs are formed 
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FIG. 2: (Color online) Total number of vortices v in the system and the sample average value \p\ of the OP as a function of 
time for 4>/4>o = 20 and a — 0.6 (a), a = 0.8(b) and a = 1.8 (c). The insets represent snapshots the amplitude of OP at t = 700 
(a) and t = 105 (b). For the quenched case (c), the snapshot displays the local vorticity at t = 20. 



without any inhomogeneity in the film. At higher fluxes 
kinematic VaV pairs are randomly created on the line 
like in the case of a "ID quench". In the x direction, 
the dynamics remains very ordered with values of the 
standard deviation of the position of the vortices V 5x 2 
approaching 0.5£. 

With the further increase of a the number of PS lines 
increases and the kinetics becomes more stochastic be- 
cause of the interaction of different PS lines. As a result, 
straight lines are replaced by vortex rivers which become 
broader and have finite curvature (Fig. 2(b)). The vortex 
rivers are comparable to the vortex self-organization dis- 
cussed in 2l| under different boundary conditions. Along 
one vortex river, few vortex-antivortex pairs are propa- 
gated. The kinetics is determined by the motion of these 
pairs along the rivers and finally by their annihilation. 
Importantly, the sample average of OP never reaches 
zero, contrarily to the case of the large field a > a C 2- The 
total number of vortices in the beginning of the process 
is larger than (Fig. 2(b)) which is also an indication 
of the growing importance of chaotic behavior in the dy- 
namics. The values of V 5x 2 are also strongly enhanced, 
reaching 2irR/3. The velocity of vortices along the rivers 
becomes smaller which is seen from the time dependence 
of the vortex number (Fig. 2(d)). Nevertheless the ve- 
locity is still high compared to the case when the OP 
has recovered to its equilibrium value. The last regime 
a > a C 2 is presented in Fig. 2(c). Here the quench con- 
dition is satisfied and the OP decreases uniformly until 
it reaches zero (Fig. 2(c)). As a result, the new phase 
starts to grow uncorrelated and the vortices are created 
randomly. The number of vortices is substantially larger 
than 4^ . Most of these vortices recombine rapidly. The 
remaining vortices move slowly through the sample prop- 
agating the 27T phase jump. The random dispersion of 
these vortices is a fingerprint of the KZ mechanism. In- 
deed, V Sx 2 reaches now 2ttR/2, which means that vortex 
distribution is completely random. Another characteris- 
tic of the KZ scenario is that the vortices are created 



while the order parameter is very close to zero and not 
during the fast growth like in the previous cases as one 
can see by comparing Fig. 2(a) and (b) with Fig. 2(c). 
It is important to notice that the total net vorticity is 
strictly equal to zero at any time in the case of periodic 
boundary condition in the z direction. 

For vacuum boundary conditions flj| the kinetics is 
very similar. When a c \ < a < a C 2 and ^ < 10 one or 
more lines with reduced OP are formed. The difference 
is that the PS lines here have finite curvature, because 
they start to grow from the edges of the film and finally 
connect each other. Further increase of the flux, keeping 
a constant leads to the formation of flux rivers. The most 
important difference is that not all "rivers" necessarily 
connect two edges of the film. As a result some of them 
ended in the middle of the film, leading to the relatively 
small vorticity. These remaining vortices and antivortices 
propagate slowly to the edges of the film and kinetic is 
determined by the slow vortex motion. The dynamics 
when a > a C 2 is governed by KZ mechanism, as in the 
previous case but the total net vorticity may be finite. 

In the case of an inhomogeneous superconductor, the 
effective coherence length is now z-dependent £(z) = 
£ / ^b(z). Therefore only the middle part of the ring may 
be unstable while the other parts of the film remain in 
the metastable state. The introduction of z-dependence 
of the parameters in Eq.([T]) is designed to enhance the 
transverse vortex dynamics and allows to demonstrate 
different mechanisms of vortex formation. As expected, 
the PS dynamics starts first in the region with stronger 
current and is characterized by a and 

In the region a c \ < a < a C 2 and small flux -4^ < 10 the 
initial stage of the kinetics is similar to kinetics in the 
homogeneous film. The VaV pairs are not well defined. 
However, when kinematic VaV pairs approach the low 
current regions, they become well defined and are slowing 
down (fig. 3 (a)). Therefore vortices are stabilized near 
the line where the tangential velocity has discontinuity. 
These vortices represent another case of KH instability in 
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FIG. 3: (Color online) Total number of vortices v in the system and the sample average value \p\ of the OP as a function of 
time for al — 0.4, <j)/<j)o = 5 and b — 4 (a); 4>/4>o = 20 and b = 2.25 (b); 4>/4>o = 50 and b — 12.25 (c). The insets represent 
snapshots of the amplitude of the OP at t = 120 (a) ; t = 790 and t = 830 (b). For the quasi quenched case (c), the snapshot 
displays the local vorticity at t = 112. 



superconductors. This instability leads to the formation 
of well defined vortices and governs the kinetics of the PS. 
To the best of our knowledge this is the only instability 
which allows vortex production in the low flux limit. 

When the flux through the ring is large > 10, the 
initial fast dynamics is similar to the dynamics in the 
homogeneous case until vortices reach the low current 
regions. Then they become slow and well defined. As it 
is seen in Fig. 3(b), the vortices propagate one after an- 
other to the film edge, demonstrating the vortex-vortex 
attraction even in the case when the OP has already re- 
covered. 

The further increase of a > a C 2 leads to the quench in 
the middle part of the film (Fig.l). During the quench 
many KZ vortices are created. Most of them are anni- 
hilated on a very short time scale. The rest reaches the 
line separating the region with different currents. The 
vortices almost stop near this line. The further dynam- 
ics is determined by the diffusion of these vortices to the 
film edges. When a is large enough, the KH vortices be- 
come well defined before the recovery of the OP in the 
middle part of the film and therefore the inhomoge neity 
suppresses the KZ mechanism in agreement with Ref. [14j , 
making kinetics less stochastic. 

Experimentally, observing such dynamics of vortices 
might be a real challenge because the short characteris- 
tic times does not allow the use of instruments with suf- 
ficient space resolution. However, recent works llj, \22 \ 
showed that freezing the dynamics can characterize both 
KZ and vortex river scenarios. Another idea is to use 
time resolved femtosecond optical spectroscopy as pro- 
posed in the Ref. (23| . As it is shown in the Refs.0, 0] 
the role of heating is not important for the proposed ge- 
ometry of the film. 

We have considered the kinetics of the flux penetra- 
tion to the 2D ring. We found out that for small values 
of the external field a, the kinetics is deterministic and 
essentially ID. Increasing the flux <f> creates kinematic 
vortices and even leads to a ID quench along the PS 



line which is a first step towards stochastic behavior (see 
Ref. [HI)- Further increase of a leads to the formation of 
vortex rivers, and ultimately to the quench of the sample 
leading to the stochastic dynamics of KZ vortices. The 
dynamics in the inhomogeneous film demonstrates that 
the VaV pairs are the topological analog of the PS mech- 
anism in 2D but this analogy is not as straightforward as 
is often believed. Finally, our calculations for a partially 
quenched film indicate that KH vortices at the interface 
are strongly predominant. 
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PHASE DIAGRAM OF THE DIFFERENT 
DYNAMICS 

We observe different dynamics in the homogeneous 
case, depending on the vaiue of the magnetic vector po- 
tential a and the number of flux quanta 4>/4>o penetrating 
the cylinder. In order to plot a phase diagram and dis- 
play the different types of dynamics, we needed to choose 
between multiple criteria. We found out that the stochas- 
ticity is a good approach to this problem and chose the 
standard deviation of the x coordinate of vortices to de- 
fine regions where the dynamics are comparable. The 
normalized standard deviation ^Jj^~ is shown in Fig.l (R 
is the radius of the cylinder). The different regions where 
the standard deviation is similar are separated by lines 
of constant standard derivation. 

• Part 1 of the phase diagram corresponds to the 
phase slip (PS) line solution. No vortices are 

present and we define there v _fg ? = because the 
dynamics is very ordered. 



ttR 



Part 2 (0 < 



ttR. 



< 0.23) corresponds to kinematic 



vortices traveling on a single line. 



/ 5x 2 



• Part 3 (0.23 < — - 
vortex rivers : stochasticity is increased by the the 



< 0.48) is assigned to multiple 



increasing number of rivers. 



Part 4 (0.55 < ^p) corresponds to Kibble-Zurek 
(KZ) type dynamics which displays the highest pos- 



sible stochasticity. 



Part 5 (0.48 < 



/ Sx 2 
nR 



< 0.55) represents the region 



where the standard deviation fails to clearly distin- 
guish between KZ and vortex rivers. In this region, 
the area with low a corresponds to vortex rivers and 
should belong to the part 3. Indeed, when the flux 
is increased at low a, the number of vortex river in- 
creases and the vortices often travel from one river 
to the other : their position becomes stochastic. 
On the opposite, at low flux, the area with high a 
should belong to part 4. 



HEATING 

The estimate of the heating of the cylinder may be 
done using Ohm's law. We consider that a part of the 




FIG. 1: (Color online) Phase diagram of the possible evolution 

drawn from the normalized standard deviation "^g 2 f the x 

coordinate of the vortices. We chose to define v^g 2 — when 
no vortices are present (Part 1). Part 1 corresponds to the PS 
line, part 2 the kinematic vortices, part 3 to vortex rivers and 
part 4 to KZ type dynamics. The part 5 represents a region 
where the standard deviation is not sufficient to distinguish 
between the dynamics. The standard deviation is plotted as 
a function of the vector potential a and of the number of flux 
quanta (a) for u = 10. a c i is the critical value under which 
dynamics are very improbable and a C 2 is the critical value 
over which we calculated KZ dynamics would occur. 



cylinder of area 7r£ 2 for the vortex or 2ttR£, for the PS 
line, behaves like a normal conductor of conductivity a n 
until its annihilation and sustains a current J s (in cgs 
units). The energy dissipated per unit of volume is, as in 



.1] 



J 2 



[<f> a(l - a 2 )]'' 

167T^A e 2 ff 



(1) 



This energy is less than a quarter of the condensation 

H 2 

energy le C2 -> ■ In the case of PS lines and kinematic vor- 
tices, when dynamics are fast and the areas behaving as 
normal conductors are small, the heat will be dissipated 
along the sample and through the contacts. It will not 
impact the dynamics. However, in the case of multiple 
vortex rivers or in the quenched KZ scenario, the total 
heating will be much larger but can be tackled by modern 
experimental techniques. Indeed, modern cooling meth- 
ods are fast enough [2|. 
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